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The characteristics of a meaningful effective field theory (EFT) analysis are dis- 
cussed and compared with traditional approaches to NN scattering. A key feature 
of an EFT treatment is a systematic expansion in powers of momentum, which 
is demonstrated using an error analysis introduced by Lepage. A clear graphical 
determination of the radius of convergence for the momentum expansion is also 
obtained. I use these techniques to compare cutoff rcgularization, two forms of 
dimensional regularization, and the dibaryon approach, using a simple model for 
illustration. The naturalness of the parameters and predictions for other observ- 
ables are also shown. 



1 Introduction 

An effective fieid theory (EFT) description of nucleon-nucleon (NN) scatter- 
ing is an important step on the road to a consistent description of nuclear 
matter. It also provides the techniques to predict other observables such as 
bound state expectation values. However_lhe interaction of two heavy nucle- 
ons requires a nonperturbative treatment Era which has lead to disa; 
about the nature and limitations of an EFT expansion in this case 
Regularization is needed to handle divergences that arise, but the results are 
said to depend on the regularization scheme used and the size of the scattering 
length involved. More generally, it is claimed that the behavior and predic- 
tive power expected from a true effective field theory is aat exhibited by every 
regularization method when applied nonperturbatively.tau 

Therefore, R. J. Furnstahl and I set out to determine the important fea- 
tures of a nonperturbative EFT^|. This can be done most clearly by making a 
side-by-side comparison of the regularization schemes listed in Table |l| using 
the error analysis advocated by Lepage. □ This comparison illustrates which 
schemes behave like a true EFT, EUl as will be seen below. In addition, the 
differences seen from a comparison with successful phenomenological models 
for NN scattering such as the Reid potentialEj help determine the importance 
of a reformulation in terms of an EFT. 
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Table 1: Regularization schemes and their abbreviations used throughout this paper. 



name regularization scheme 




CR[G] Cutoff regularization with gaussian weighting in the potentiall 
DR[MS] Dimensional regularization with modified minimal subtraction. 
DR[PDS] Dimensional regularization with power divergence subtraction!! 
dibaryon DR[MS] but with an additional degree of freedom associate) 
with a low-energy bound or nearly bound state Lni^. 
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2 Comparison Techniques 

1 first describe how to analyze and compare the different regularization schemes 
in a way that will allow for a determination of which schemes have the system- 
atics and predictability desired of an effective field theory. Concentrating on 
S-w&ve scattering, I will focus only on the short-distance EFT. This general 
5-wave potential can be written as 



V(p, p ) = — c - OA9 + e — + . . . , (1) 




with a scale A s introduced to make the coefficients c, d, e, . . ., dimensionless. 
There is also a scale A implicit in this effective potential that signifies the point 
above which new degrees of freedomj-not accounted for in the effective theory 
begin to contribute to the physics. eI The factor of l/4n picked up by each 
additional term in the Born series requires a An to be factored out in order 
to render these coefficients natural, q i.e. of order unity, as discussed further 
below. 

The coefficients are determined order-by-order from matching to the avail- 
able data. To fit the first n constants in the potential requires at least n points 
of data. The potential is not a measurable quantity, so instead a scattering 
observable such as the phase shift S(p) is used to determine these constants. 

Inserting this potential into the Lippmann-Schwinger equation and spec- 
ifying a regularization scheme, the resulting amplitude determines the phase 
shift by one of the following two equivalent relations (T = —A) 

T(p) = -^e iS ^smS(p) , "Sr^ =P™t6(p) - ip . (2) 

For example, dimensional regularization with modified minimal subtraction 
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(DR[MS]) gives 



T V & p) (3) 

1 BR[MS] ~ iM p , * W 
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The imaginary part of the denominator is linked to the numerator by unitarity. 
In general, however, the real part of the denominator can be more complicated 
than in the DR[MS] case and consequently dictates the behavior of the scheme 
near a bound state. This can be seen by calculating the same quantity for 
dimensional regularization with power divergence subtraction (DR[PDS]) 

_ V(p,p) u , 
Jdr[pds] jj ( 4 ) 

47T 

which reduces to the DR[MS] result for \l = 0. However, the appearance of 
this additional mass scale is important for a good radius of convergence as will 
be seen in the next section. 

Cutoff regularization is another method for dealing with divergences in 
which momenta greater than the scale of new physics A are explicitly sup- 
pressed. This can be implemented by adding a gaussian weight exp(— (p 2 + 
p /2 )/2A 2 ) to the potential in Eq. (|l]) and will be referred to as CR[G]. Although 
a physically intuitive method, numerical techniques are required to solve this 
type of regularization. However, one can see the effect of the cutoff by taking 
a geometric sum of the first two Born terms of the CR[G] method. This just 
ends up giving Eq. (0) with /i = 2A c /n. It will be shown below that DR[MS] 
has very different results from both DR[PDS] and CR[G], which can be traced 
to this difference in the real parts of the T-matrix denominators. 

In the left plot of Fig. [tj I show the results from using CR[G] to fit one, 
two, and three constants in the potential, as compared to the exact S-wave 
phase-shift for the delta-shell potential, EJ which models the underlying physics 
and is discussed in detail in the next section. A first glance at the plot shows 
the approximation to the phase shift improves as more constants in the po- 
tential are fit. However, a second look shows that it is very difficult to gather 
any quantitative information from the plot. At what point the curve deviates 
enough to be considered inaccurate is not clear, and the radius of convergence 
of the EFT expansion is completely obscure. Furthermore, a simple calculation 
shows that every termJn the momentum expansion of the phase shift contains 
the scattering length. Q This means that the error in the phase shift could be 
numerically sensitive to a large scattering length and contaminate the power 
counting. 
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Figure 1: The phase shift <5(p) (left) and the error in pcot<5(p) (right), each plotted as a 
function of p for the delta-shell potential with a weakly bound state, as discussed in Section 3. 
The solid line is the exact result and the dashed lines show the CR[G] fit for one, two, and 
three constants. 



This issue is avoided with no loss of generality by plotting the error of 
pcot 5{p) instead: 

A[p COt S(p)} = |p COt tf eff (p) -pCOt£ t rue(p)| • (5) 

If the effective field theory follows proper power counting, then this error should 
improve by two powers of momentum as each additional coefficient in the po- 
tential is fixed in Eq. (|l|). Since a large scattering length is synonymous with 
a near bound state (or pole) in the amplitude, Eq. (||) shows that this pole 
is cleanly mapped only to the first constant in a momentum expansion of 
p cot S(p). It is known from conventional scattering theory that this combina- 
tion has a well defined expansion in p 2 for short-range potentials known as the 
"effective range expansion" 

pcotS(p) = -— + \r e p 2 + v 2 p k + . ■ . , (6) 
a s 2 

which defines the scattering length a s and effective range r e . When long-range 
potentials are included, Eq. ra) is only valid at low momentum or is even 
inapplicable. In this case, one must define a modified effective range expan- 
sion. E3 Since the effective theory contains the same long-distance physics as 
the true underlying theory, Eq. (|^) can be modified to have a clean momentum 
expansion. For short range potentials as considered here, Eq. (0) is sufficient. 
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Plotting the error in p cot 5{p) as a function of p on a log-log plot, a straight 
line is_expected with slope given by the dominant (lowest) power of p/A in the 
error. □ As more constants are included, the slope in this error should increase, 
signifying the removal of higher powers of p/A. The second plot of Fig. [j] 
clearly demonstrates the order-by-order improvement in the amplitude as more 
constants are added to the effective potential. With one constant the slope of 
the error is two, i.e., 0(p 2 / A 2 ), and increases by two with each additional 
constant. 



3 Illustration with Delta-Shell Potential 

In this section, I illustrate the points made above by using EFT techniques with 
the different regularization methods to systematically describe the "unknown" 
short-range physics of a specific example. I could use NN scattering data, 
but it is more convenient to use an exactly solvable potential to serve as data 
in order to have a clean understanding of what features are important. The 
delta-shell potential has been used in the past to simulate the large scattering 
length found in NN scattering. EJ-Kaplan used this potential to illustrate the 
benefits of the dibaryon approach. He found upon considering NN scattering 
that the inclusion of long-distance pion physics dkLnot change the conclusions. 
This agrees with the experience of other authors Bel that the addition of pions 
as long-range interactions does not affect the properties of the short-range 
expansion. The delta-shell potential is therefore a sufficient model for the 
purposes here. 

The potential can be written in terms of the nucleon mass M, the coupling 
<7, and the range of the potential r . This short-range potential represents the 
new physics of the underlying model theory, and so r — 1/ A is taken below, 

VtMr) = -g^s(r-]-). (7) 



M V A , 

It has exactly one bound state for g > 1 and no bound states for g < 1. 
Scattering with p > A probes the details of the potential, so one would expect 
A to be the radius of convergence of a well-tuned EFT. The scattering length 
becomes very large for g near 1, whereas the effective range (and the rest of 
the terms in the effective range expansion) are of natural size for all g: 

„ _ 9 1 2( 5 + l) 1 



5-1 A 3g A 

The scattering length in the 1 S'o channel of iViV-scattering can be modeled by 
choosing (g, A) = (0.99, m p ). This potential Eq. (0) with different g's will be 
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Figure 2: The error in p cot 5(p) plotted as a function of p/A for a small scattering length 
without a bound state g = — 10 and for a large scattering length with a bound state g = 1.01. 



the "laboratory" from which the different regularization schemes are compared. 
I will also briefly discuss results for actual measured data and the inclusion of 
pions below. 

First, some intuition will be gained by graphically reproducing Kaplan's 
result that DR[MS] has a small radius of convergence if the scattering length 
is large. From Eq. (ph, it can be seen that choosing g — 0.99 gives a scattering 
length one hundred times larger than choosing g = —10. At the same time, the 
effect from the presence of a bound state is investigated by taking g — 1.01. 
The effective potential is given by Eq. ([!]) with the mass scale A s associated 
with the inverse delta-shell radius and the prescription of using DR[MS] on all 
divergent integrals. 

The momentum expansion for DRPMS] can be shown analytically l3'Ll to 
be p 2 a s r e /2. This implies using Eq. (g) that the radius of convergence for 
g = 1.01 and 0.99 should be roughly 1/10 that of the g = — 10 case. Fixing the 
constants in the potential Eq. (|l|) by matching p cot 8(p) produces the results 
in Fig. ||. The dashed lines show the DR[MS] results for one, two, and three 
constants respectively. Indeed all three lines converge to p/A ~ 1 for g = — 10 
and p/A ~ 0.1 for g = 1.01. The results for g — 0.99 fall on top of the g = 1.01 
results and are therefore not shown. This implies the presence of a bound state 
does not matter, but the size of the scattering length does. It would be difficult 
to draw these conclusions had only the phase shift itself been plotted (the left 
plot of Fig. |). Of course, these results can also be shown analytically for this 
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Table 2: Effective potential for g = — 10 (small scattering length) to three different orders 
in the error 0(p n / A n ) for different regularization schemes. Aa = 1 for CR[G] and fi = A for 
DR[PDS]. 



error 
order 


c 


DR[MS] 

d e 


c 


DR[PDS] 

d e 


c 


CR[G] 
d e 


2 


0.758 




8.33 




1.49 




4 


0.758 


-0.206 


8.33 


-25.0 


2.67 


-1.04 


6 


0.758 


-0.206 0.0672 


8.33 


-25.0 38.2 


2.67 


-1.58 3.75 



Table 3: Effective potential for g = 1.01 (large scattering length) to three different orders in 
the error C(p n /A n ) for different regularization schemes. Aa = 1 for CR[G] and fi = A for 
DR[PDS], 



error 
order 


DR[MS] 

c d e 


DR[PDS] 

c d e 


CR[G] 

c d e 


2 
4 
6 


84.2 

84.2 -5640 

84.2 -5640 378000 


-0.842 

-0.842 -0.564 
-0.842 -0.564 -0.152 


-1.42 

-0.946 0.842 
-0.937 0.618 -0.181 



simple model, but the analysis applies much more generally, when part or all 
of the calculation is done numerically. 

The constants c, d, and e are given in Tables ^ and 0. Their values are 
required to at least 8 digits to produce the accuracy of Fig. |[ The correlation 
between the naturalness of the constants and the radius of convergence is 
apparent. The constants for the g = 1.01 case are extremely large, reflecting 
the breakdown of the effective field theory for DR[MS] much below the expected 
scale A. 

One way to fix this behavior is to introduce the dibaryon.Em This takes 
the large scattering length into account by explicitly introducing a low-energy 
s-channel degree of freedom into the effective lagrangian. For g > or g < — 1, 
the potential can be written as 

v , f\ V 2 r < 277 2 SnM + g 3A 2 1 - g 

(9) 

with E — p 2 /M always kept on-shell. Since there seem to be three constants 
fit in Eq. (|), one might think the dibaryon will have an error of C£p 6 /A 6 ). 
However, the dibaryon amplitude is only matched to second order u in the 
momentum when deriving the relations in Eq. @ and indeed shows an error 
of 0(p 4 /A 4 ) for the two values of g in Fig. || (plotted as the dot-dashed line). 
The slope and magnitude of the error do not depend on the scattering length, 
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as expected. 

I next repeat the calculations using the cutoff regularization method CR[G] 
with A c = A. There are various ways to numerically solve the Schrodinger 
equation with a cutoff, but the following procedure is particularly efficient and 
numerically robust enough to attain the accuracy required here. First, the 
variable phase method is used to solve for the phase shift. This is a differential 
equation, 

S'(r) =-— V(r) sm 2 (pr + 8(r)) with 8(0) = , (10) 
V 

which expresses the change in the phase shift as the potential is built up from 
zero at r = to its full value at r — oo. The boundary condition ensures that 
the full phase shift given by S(oo) is zero in the absence of a potential and 
defines the otherwise ambiguous multiple of tt in the phase shift. The routine 
ODE from package ODEE3 was used to solve the differential equation. 

Then I use a general method for taking the theoretical error of the EFT 
into account when matching the constants in the effective potential. After 
evaluating the combination p cot S e g(p) in the effective theory, it is subtracted 
from the true result (either data or an exact solution to a model problem) and 
the difference, 

2 4 

ApcotS(p) = a + /?^+ 7 | I + ... , (11) 

is fit to a polynomial in p 2 /A 2 to as high an order as possible. To obtain the 
accuracy shown in the plots, a weighted polynomial fit of Ap cot S(p) up to 
p/A = 0.1 was required. 

Using a spread of momentum near zero, the polynomial fit should be 
weighted with both the expected theoretical error in momentum and any ad- 
ditional experimental noise. The resulting coefficients a, /?, 7, . . . are then min- 
imized with respect to variations in the effective potential constants c, d, e, . . . 
using an optimization code. In practice, this method is more robust and nu- 
merically stable than matching the values of p cot 8 Lp) at discrete points to 
fix the constants. This allows the analysis of LepageQ to be extended beyond 
second order. Also note that such a procedure is needed when matching to 
data even when the EFT observables can be calculated analytically. 

The number of coefficients that can be minimized is given by the number 
of constants retained in the effective potential. Liised DPOLFT from package 
SLATEC0 to find the polynomial fit and MINF,0 which is based on the Powell 
method, to carry out the minimization. Normal accuracies in minimization 
using double precision numbers with this method are 10~ 12 or better. 
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Figure 3: The error in p cot 5(p) plotted as a function of p/A for a small scattering length 
without a bound state g = — 10 and for a large scattering length with a bound state g = 1.01. 

The solid lines in Fig. || show that CR[G] does work regardless of the 
scattering length. In fact, with A c ~ A, the result is just as good as the 
dibaryon for the same number of constants. The values for the constants are 
given in Tables || and ||, showing they are all natural for both g's (although 
the third constant is somewhat small for g = 1.01). 

Note that as more constants are fixed, the lower-order constants are mod- 
ified. This occurs because even after truncating the potential for the cutoff 
regularization to a given order, it still contains all orders in p 2 from the gaus- 
sian factor. The nonperturbative solution of the Lippmann-Schwinger equa- 
tion therefore can generate terms of any order in p 2 /A 2 . The amplitude itself 
is matched to the true result order-by-order in the momentum so the power 
counting of the potential is destroyed. This consequence of the cutoff regular- 
ization is not necessarily relevant since the potential is not an observable. It 
is interesting to note, however, that these modifications are relatively small. 
This is also true when adding the long distance physics of the pion in fitting 
to actual NN data as shown in the next section. 

In summary, Fig. ^ shows that all regularization schemes considered so far 
produce useful effective field theories for a small scattering length, but DR[MS] 
fails for large scattering length. A failure of the power counting in powers of 
p/A is reflected in unnatural constants in the potential. _ 

I now turn to the most recently proposed regularization scheme,li3 DR[PDS]. 
An additional prescription compared to the DR[MS] case is an expansion of 
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observables to the same order in p 2 as the potential Eq. (|l|). If this is not done, 
the results are fi dependent, with fj, = reproducing DR MS] and \i larger than 
the nucleon mass approaching the CR[G] result in Fig. |. 

When only dealing with a short-range potential as discussed so far, the 
DR[PDS] prescription reproduces the effective range expansion Eq. @ by con- 
struction. The DR[PDS] results in Fig. | are therefore fj, independent, but the 
constants are still /j, dependent. I take ^ = A to compare with CR[G]. This 
produces natural constants for g = 1.01 but somewhat unnatural ones for 
g = —10 as seen in Tabic |[ This is an accidental consequence of the momen- 
tum expansion being in powers of p 2 r e /2(/i — l/a s ), so that for g = —10 and 
/x = A, Eq. (||) shows the denominator is nearly zero. This implies that the 
scale (x is not functionally equivalent to the cutoff 1/a in CR[G] since it does 
not always signal the onset of new physics at the scale A. 

For both large and small scattering length, DR[PDS] does quite well, with 
a radius of convergence p/A ~ 1. The CR[G] result is better for one constant 
since the cutoff generates an effective range r e close to the true result. Overall, 
DR[PDS] is a convenient method to produce reasonable analytical results, and 
depending on the problem at hand either CR[G] or DR[PDS] may be suitable. 
One should note, howevex, that only DR[PDS] provides a strict diagram by 
diagram power counting. E3 ._. 

Finally, I mention the Reid potential^ which consists of a sum of Yukawa 
interactions, with fixed masses chosen as integer multiples of the pion mass, 
and coefficients that can be varied to fit the data. So for S'-waves, 

wp, po = + q2+ (3 m ,)2 + q2+ (L.)2 + • • • ( i2 ) 

= opF c 2 (5m T ) 2 + e 3 (3m 7r ) 2 + (c 2 + c 3 )q 2 
[q 2 + (3m.) 2 ][q 2 + (5m.) 2 ] 

Combining the Yukawa terms in the second equation generates a numerator 
that looks like an EFT and a denominator that suppresses large momentum, 
similar to the cutoff regularization. However, since each term has a different 
mass, a clear separation scale A is not identified. There is also a one-pion 
exchange (OPE) term which is not important for the short-range physics dis- 
cussed here. 

The original Reid analysis used a global fit and only one adjustable param- 
eter in each channel, for a result with approximately the same error (roughly 
a few percent of the data) at all momenta. However, the EFT fitting proce- 
dure can be applied instead. Using Yukawa masses comparable to A should 
give similar results to CR[G]. Indeed, if the OPE contribution is ignored and 
a low-momentum fit is done to the constant C2, the error plot is similar to 
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the CR[G] result with one constant (Fig. ^). Adding a second short-range 
Yukawa does as well as CR[G] with two constants, since the Yukawas play off 
each other to allow the next order error in q 2 to be removed. This interplay 
becomes increasingly complex at higher orders. Furthermore, the additional 
mass scales obscure (or smear out) the role of A as a scale that separates the 
known from the unknown physics in effective field theories. Traditional NN 
potentials such as Reid are well suited for global fits, but systematic predictions 
with controlled error estimates are more properly analyzed using an effective 
field theory. 

4 Other Observables 

I now turn to an investigation of the binding energy. If the EFT is truly repro- 
ducing the S-matrix of the underlying theory order-by-order in a momentum 
expansion, it should reproduce the binding energies and other observables to 
the same order of accuracy as the phase shifts. I therefore use the binding en- 
ergy prediction as a consistency check for the candidate effective field theories 
found so far. 

The potentials have already been fit to a given order by the scattering 
phase shifts above, and I use these potentials without adjustment to solve for 
the binding energy. This can be done analytically for the DR schemes by 
finding the poles in the scattering amplitude. The exact binding energy -Ebind 
for the delta-shell potential is given by solving the equation 

1 1-e- 2 " VME bind 

- 9 = ^r> v (13) 

There is only one bound state to predict in the delta-shell potential, and if it is 
shallow enough, even the effective range expansion with the values of a s and r e 
can determine its value. A better test is to increase g until -Ebind is large and 
on the order of A, and use the EFT to determine the accuracy of the binding 
energy prediction as a function of this variation. If a true radius of convergence 
is present, the effective field theory should break down for E'bind/A ~ 1. The 
binding energy is 6.48 x 10 -2 MeV for g — 1.01 but quickly increases to 812 
MeV for g = 2.5. 

The relative error in the binding energy is plotted in Fig. [|. Both CR[G] 
and DR[PDS] show the clear power counting behavior and proper radius of 
convergence expected from a true effective field theory. This gives a graphical 
verification that the errors in the binding energy really do follow power count- 
ing rules. I have checked that the same behavior is seen when plotted as a 
function of the average momentum \J (p 2 ) /A. The dibaryon result also follows 
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Figure 4: The error in the binding energy for the candidate effective field theories for repre- 
sentative values of g from 1.01 to 2.5. The star signifies the absence of real binding energies 
for DR[MS] with two constants at the values of g > 1.25 considered. 

the expected error scaling. In contrast, the deficiencies of DR[MS] regulariza- 
tion seen for the phase shifts are manifested here as binding energies that do 
not follow the EFT error scaling, improving to a lesser degree than expected. 
In addition, for values of g > 1.25 with two constants, the DR[MS] S-matrix 
shows no bound state with a real energy. 

Therefore, these results show that most regularization procedures demon- 
strate the characteristics of a systematic predictive effective field theory. The 
fit of more and more constants in the effective potential improves the predictive 
power order-by-order in the momentum expansion. The radius of convergence 
of the EFT is independent of the scattering length and is given by the scale 
where new physics enters. 

This test can be taken further by looking at moments of the bound state 
wavefunction, such as (r 2 ). Preliminary results are shown in the left plot 
of Fig. |^ for one and two constants. Going touihird order is not productive 
here because the effective operator r 2 changes □ due to the operator product 
expansion at this order in 1/A 2 . A systematic improvement of the expectation 
value is nevertheless seen for one and two constants. The radius of convergence 
does not look optimal, but this is merely an illusion since CR[G] with one 
constant does slightly better than anticipated for large i?bmd/A, as seen from 
Fig-!- 

Finally, it would be ideal to apply the present analysis to real data for 



12 




Figure 5: The left plot shows the error from the bound state moment (r 2 ) in the delta-shell 
potential. The right plot shows the result for one, rtyo, and three constants in CR[G] when 
fit to the 1 5o NN scattering data from Nijmegen. La 



NN scattering. Using the data from the Nijmegen phase shift analysis, E3 
I have taken the theoretical error in the EFT analysis with CR[G] and the 
experimental uncertainties into account in constructing the right plot of Fig. ^. 
Since two-pion physics was not included in the analysis, the cutoff was taken 
near the 2to w threshold (A = 250 MeV). The one and two constant results 
show the systematic improvement and nice radius of convergence expected 
from above. However, the uncertainty in the long distance physics blurs the 
three constant result given by the range between the two dotted lines. This 
must be looked at in more detail to achieve a better precision fit. 



5 Conclusion 

I have presented resultsS showing how the error analysis suggested by Lepagei 
can be applied to compare cutoff regularization, two forms of dimensional 
regularization, and the dibaryon approach in the context of nonperturbative, 
nonrelativistic effective field theories. This analysis focuses on a key signature 
of EFT behavior: the systematic scaling of errors with momentum or energy. 

New numerical procedures allow for an analysis to third order and be- 
yond in the EFT expansion, which is necessary to obtain a clear graphical 
determination of the radius of convergence for a given observable. Such an 
analysis is required for a systematic fit to data regardless of the regularization 
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scheme. 

It is found that all of the regularization methods except for dimensional 
regularization with modified minimal subtraction are consistent with basic 
features expected from a useful effective field theory: 

• Each additional order in the potential leads to a systematic improvement 
in the amplitude. 

• The radius of convergence for this improvement, when optimized, is dic- 
tated by the scale of new physics. 

• Other observables are predicted with the same accuracy as the amplitude 
at each successive improvement. 

These results are consistent with the analysis of van Kolck that, with proper re- 
summations, any effective field theory for short-range interactions is equivalent 
to an effective range expansion to the same order. Ilia 

The CR[G] and DR[PDS] regularization schemes are each suitable for de- 
veloping effective field theories of many-nucleon systems. In future work both 
schemes will be used in extending this fitting procedure and error analysis to 
nuclear matter. 
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